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Long-term monitoring of LS I +61°303 with INTEGRAL 
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ABSTRACT 

LS I +61°303 is one of the few high-mass X-ray binaries that have been recently- 
observed at TeV and GeV energies. Here we investigate the hard X-ray spectral and 
timing properties of this source using the IBIS/ISGRI instrument on-board the INTE- 
GRAL satellite. We carry out a systematic analysis based on all available INTEGRAL 
data since December 28, 2002 up to April 30, 2008. The total exposure time analyzed 
amounts to 2.1 Ms, hence more than doubling the previous reported sample. The 
source is best detected in the 18-60 keV band, with a significance level of 12. Oct. The 
hard X-ray data are best fit with a simple power law with a photon index of ~ 1.7±0.2. 
We detect a periodical signal at 27 ± 4 days, matching the orbital period of 26.496 
days previously reported at other wavelengths. The hard X-rays orbital lightcurve is 
obtained and compared with those derived at other frequencies. 
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1 INTRODUCTION 

LS I +61°303 is among the few binaries observed to date 
to emit high energy 7-rays. It hosts a BO Ve main-sequence 
star (Hutchings & Crampton 1981; Casares et al. 2005), or- 
bited by a compact object of yet unknown nature, with a 
period of 26.4960±0.0026 days (Gregory 2002). The system 
is located at a distance of about 2 kpc. The orbital zero 
phase is taken at JD 2443366.775 (Gregory & Taylor 1978), 
and the phase values of 0.23±0.02 (Casares et al. 2005), 
-0.275 (Aragona et al. 2009), and 0.30±0.01 (Grundstrom 
et al. 2007) represent the periastron passage uncertainty. 

LS I +61°303 has been recently detected as a periodic 
7-ray source by the Major Atmospheric Imaging Cerenkov 
telescope (MAGIC, Albert et al. 2006) and The Very En- 
ergetic Radiation Imaging Telescope Array System (VERI- 
TAS, Acciari et al. 2009). A distinctive orbital modulation 
of the VHE 7-ray emission was seen, which was found to be 
anti-correlated with that observed at GeV energies by the 
Fermi gamma-ray satellite (Abdo et al. 2009). In the last 
few years, there has been a burst of activity trying to un- 
derstand the nature of this source, whether it is composed 
by a pulsar or a black hole system, and which are the mech- 
anisms that lead to their multiwavelength behavior (e.g., 
Dubus 2006, Gupta & Bottcher 2006, Romero et al. 2007, 
Sierpowska-Bartosik & Torres 2009). 

The current observations on LS I +61°303 are rather 
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sporadic, especially at high energies, which prevent from 
monitoring the phase evolution on time scale of years. The 
limited observations at soft X-rays conducted by XMM- 
Newton (Neronov & Chernyakova 2006; Sidoli et al. 2006), 
Chandra (Paredes et al. 2007, Reaet al. 2010), ASCA (Leahy 
et al. 1997), ROSAT (Goldoni & Mereghetti 1995; Taylor 
et al. 1996) and Einstein (Bignami et al. 1981) were too 
short to cover an entire orbit. An exception to this were the 
72XTE/PCA observations. 7?XTE observed LS I +61°303 
in 1996 March 1 - 30, and found that the 2-10 keV flux 
peaks at the orbital phase ~ 0.45 - 0.6. The detection and 
non-detection by HEXTE suggest a similar trend in orbital 
modulation at 15-150 keV (Harrison et al. 2000). Smith et al. 
(2009) adopted about half year of 7JXTE/PCA observations 
covering the time period from August 28, 2007 to February 
2, 2008 (MJD 54340-54498) and reported a possible two- 
peak orbital lightcurve at 2-10 keV band and several very 
short flares (which due to the large FoV of the PCA instru- 
ment, could still come from other sources in the field). The 
orbital phase derived at soft X-rays (0.3-10 keV) peaks at 
the phase 0.65, close to the apastron (Esposito et al. 2007). 
This is consistent with the contemporary observations at 
above hundred GeV by MAGIC (Albert et al. 2009,2009b). 

At hard X-rays, although the source was monitored by 
INTEGRAL (Chernyakova, Neronov & Walter 2006), data 
analyzed till now was not enough to explore the appear- 
ance of the orbital period in the power spectrum nor to 
construct a detailed lightcurve. Here we investigate the hard 
X-ray spectral and timing properties of this source using the 
IBIS/ISGRI instrument on-board the INTEGRAL and carry 
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out a systematic analysis based on all available INTEGRAL 
data since December 28, 2002 up to April 30, 2008. The total 
exposure time analyzed amounts to 2.1 Ms thus enhancing 
the previously published report (Chernyakova, Neronov & 
Walter 2006) by taking into account more than 3 years of 
additional data, more than doubling the previous sample. 



2 OBSERVATIONS AND DATA ANALYSIS 

INTEGRAL (Winkler et al. 2003) is a 15 keV - 10 MeV 7- 
ray mission. Its main instruments are the Imager on Board 
the INTEGRAL Satellite (IBIS, 15 keV - 10 MeV; Uber- 
tini et al. 2003) and the SPectrometer onboard INTEGRAL 
(SPI, 20 keV - 8 MeV; Vedrenne et al. 2003). These in- 
struments are supplemented by the Joint European X-ray 
Monitor (JEM-X, 3-35 keV) (Lund et al. 2003) and the Op- 
tical Monitor Camera (OMC, V, 500 - 600 nm) (Mas-Hesse 
et al. 2003). At the lower energies (15 keV - 1 MeV), the 
CdTe array ISGRI (Lebrun et al. 2003) of IBIS has a better 
continuum sensitivity than SPI. The satellite was launched 
in October 2002 into an elliptical orbit with a period of 3 
days. Due to the coded-mask design of the detectors, the 
satellite normally operates in dithering mode, which sup- 
presses the systematic effects on spatial and temporal back- 
grounds. The INTEGRAL observations were carried out in 
the so-called individual SCience Windows (SCWs), with a 
typical time duration of about 2000 seconds each. In this 
work, only IBIS/ISGRI public data are taken into account; 
the source is too weak to be detected by JEMX and SPI. 
The available INTEGRAL observations, when LS I +61° 303 
had offset angle less than 14 degrees, comprised about 875 
SCWs, adding up to a total exposure time of ~2.1 Ms 
(covering rev. 25 - 667, MJD: 52636 - 54586). This to- 
tal exposure then enhances the previously published re- 
port (Chernyakova, Neronov & Walter 2006) by taking into 
account more than 3 years of additional data, i.e., more 
than double the previous sample. The data reduction was 
performed by using the standard Online Science Analysis 
(OSA), version 9.0. The results were obtained by running 
the pipeline from the flowchart to the image level, and the 
spectrum was derived using the mosaic images, as are ap- 
propriate for spectral analysis of faint sources. The spectra 
were fitted with XSPEC vl2.3.1 and the errors on the model 
parameters were estimated at 90% confidence level. 

In order to search for a periodic signal in the lightcurve 
data, we used the Lomb-Scargle periodogram method (Lomb 
1976; Scargle 1982) and followed the procedure described 
in Farrell et al. (2009). Power spectra were generated for 
the lightcurve using the PERIOD subroutine (Press & Ry- 
bicki 1989). We set oversampling and high-frequency fac- 
tors, which determine the period range and resolution (Press 
& Rybicki 1989), to two and four, respectively, allowing to 
search for periodic variability in the 1.1 - 1000 d range. 
The 99% white noise significance levels were estimated using 
Monte Carlo simulations (see e.g. Kong, Charles & Kuulk- 
ers 1998). The 99% red noise significance levels were esti- 
mated using the REDFIT subroutine, which can provide the 
red-noise spectrum via fitting a first-order auto-regressive 



Table 1. The flux and significance of LS I +61°303 obtained from 
the mosaic images of the combined ISGRI/IBIS observations. 



Energy 


flux 


sig. 


(keV) 


ct/s 


mCrab 


a 


20-40 


0.21±0.02 


1.60±0.16 


10.0 


40-100 


0.16±0.02 


2.10±0.28 


7.6 


100-300 


0.07±0.02 


5.00±1.25 


4.0 


18-60 


0.32±0.03 


1.60±0.13 


12.0 



process to the time series (Schulz & Mudelsee 2002, Farrell 
et al. 2009)Q 

The most significant INTEGRAL detection of 
LS I +61° 303 is derived from combining all the IS- 
GRI data, and is found at a significance level of ~ 12.0a 
in the 18-60 keV band (see Figure 1). This significance 
is consistently higher than that reported by Chernyakova, 
Neronov & Walter (2006, of 8.1a) at the 22-63 keV, based 
on 600 SCWs (covering rev. 25-288, from 2003 January to 
2005 March, and producing a corrected exposure of 273 
ks; as opposed to ours of about 578 ks). Table 1 shows the 
details on the imaging results in 4 energy bands: 20-40, 
40-100, 100-300, and 18-60 keV. The neighboring source 
QSO 0241+622 is detected at 3.4 mCrab level in 18-60 keV, 
implying a ~ 25a detection. This source is only about 1.4 
degrees away and it is usually a problem for the analysis 
of LS I +61°303 with non-imaging detectors. The energy 
spectrum of the source is produced combining all the ISGRI 
data in the 13 keV - 1.0 MeV band. Since we note that 
the sky region around LS I +61° 303 is rather crowded and 
LS I +61°303 itself is relatively faint in hard X-rays, it is 
necessary to obtain the energy spectrum from the mosaic 
images (see also Chernyakova & Neronov 2008). The energy 
spectrum is well fitted by a power law shape, with a photon 
index of 1.7+0.2 and normalization of 3.1I31 x 10~ 3 ph 
cm" 2 s" 1 keV" 1 . The flux we derive in the 20-200 keV 
band is 2.8 xl0 _11 ergs cm -2 s _1 . The reduced \ 2 f° r the- 
fit is 0.4 under 7 dof. The orbitally separated results of the 
spectral fitting can be seen in Table 2. These results show 
that our spectral fittings are mostly consistent with those 
published previously by Chernyakova, Neronov & Walter 
(2006). They also present a large improvement, especially 
for the phase interval 0.6-0.8, where the error bar in the 
photon index is largely reduced from > 1 to 0.3. 

The ISGRI lightcurve of the source is obtained in the 
18-60 keV range with a time bin following the SCWs (see 
Figure 2). A fit to this lightcurve with a constant produces 
a value of \ 2 ~ 1-2 (1058/874), indicating a deviation at a 
significance level of 4a. We then investigated whether there 
is a recognizable orbital variability in this data using timing 
analysis. The orbital lightcurve is then made by using the 
radio period and a reference time at MJD=43366.275 (Gre- 
gory 2002) at the hard X-rays, so that the absolute orbital 
phases can be directly compared with those at other frequen- 
cies (and we do this below). We find that at the hard X-rays 

1 See ftp://ftp.ncdc.noaa.gov/pub/data/paleo/softlib/redfit 
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Table 2. Orbitally separated spectra from ISGRI/IBIS observations (the results of our analysis are marked with a star) and a comparison 
with the previous work by Chernyakova et al. (2006); quoted in the three last columns. 



orbital V* (1 a) flux* (20-60 keV) effective T (1 a) flux (20-60 keV) effective 

phase (10 — 11 erg cm -2 s —1 ) exposure* (ks) (10 — 11 erg cm -2 s — 1 ) exposure (ks) 



0.4-0.6 1.9 ± 0.2 2.84± 0.54 127 1.7±0.4 3.8±0.6 50 

0.6-0.8 1.5 ±0.3 2.09± 0.52 144 3.6^1.1 3.0±1.0 23 

0.8-0.4 1A to.t L07 ± °- 36 307 1.4±0.3 2.4±0.3 200 

whole orbit 1.7 ± 0.2 1.74± 0.26 578 1.6±0.2 2.5±0.3 273 




Figure 1. The mosaic image of the LS I +61°303 sky region 
derived by combining all ISGRI data at the 18 - 60 keV band. 
The strongest source is QSO 0241+622 (at the upper part of 
the image) and the relatively faint one is LS I +61°303. The 
significance level is given in the color scale. The contours start at 
a detection significance level of 5 cr, with a step of 3 a. 
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Figure 2. The ISGRI lightcurve at 18-60 keV. The bin size is 
one SCW and the error bar is \a. 



the maxima is around phase 0.5-0.6 whereas the periastron 
of the system is at around 0.25-0.3 (Casares et al. 2005, 
Grundstrom et al. 2007, Aragona et al. 2009). The Lomb- 
Scargle power spectrum based on the ISGRI lightcurve is 
presented in Figure 3. The maximum of this power spec- 
trum is located at around 27 days with a power density of 
12.8 and an error of 4 days based on Monte Carlo simula- 
tions. This periodical signal matches the orbital period of 
26.4960 days, as reported in radio (Gregory 2002). To ob- 
tain the error in the period determination we have sampled 
the light curve that was used to produce the power spec- 
trum. The flux and flux error of each bin of this light curve 
are regarded as the mean and the deviation of a Gaussian. 
Accordingly, for each bin, a Monte Carlo flux is obtained 
to form a new light curve. In the power density of each so- 
generated lightcurves, the peak around 26.5 days is taken. 
The width of this peak distribution provides one sigma error 
of the measured orbital period 26.72 days from the real IS- 
GRI ldata. Over-plotted in Figure 3 are the white noise and 
the red noise both at 99% confidence level. The comparison 
of our INTEGRAL lightcurve with those obtained at other 
frequencies is shown in Figure 4, which is further discussed 
next. 
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Figure 3. The Lomb-Scargle power spectrum based on the ISGRI 
lightcurve (18 - 60 keV). The white noise (dashed line) and the 
red noise (solid line) at the 99% confidence level are plotted. The 
LS I +61° 303 orbital period of 26.496 days is indicated with an 
arrow. 
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Figure 4. The orbital lightcurves in a multi-frequency context, from X-rays to GeV (Abdo et al. 2009) and TeV gamma-rays (Albert et 
al. 2009). All the plots in this paper have adopted the periastron passage as 0.275 and apastron passage 0.775 as shown by the vertical 
solid lines (Aragona et al. 2009). 



3 DISCUSSION 

We have carried out a systematic analysis based on all 
available INTEGRAL data since December 28, 2002 up to 
April 30, 2008. The total amount of time analyzed was 
~2.1 Ms, thus enhancing the previously published report 
(Chernyakova, Neronov & Walter 2006) by taking into ac- 
count more than 3 years of additional data, more than dou- 
bling the previous sample. Corrected exposures after cuts 
changed from 273 to 578 ks. The difference in the amount of 
data analyzed allows for a much more detailed INTEGRAL 
lightcurve to be constructed, as can be visually seen in 
panel 3 of Figure 4 above, where we show the Chernyakova, 
Neronov & Walter (2006) 20-60 keV lightcurve in the back- 
ground, and our current results superposed. Our analysis of 
INTEGRAL /ISGRI data results in the first direct orbital 
period determination through a power spectrum analysis of 
the hard X-ray data of LS I +61° 303, although not yet truly 
significant, and in a spectrum at hard X-rays which is simi- 
lar to the previous report (Chernyakova, Neronov & Walter 
2006), but with detection significance improved to 12. Oct, 
meaning that the source did not show significant changes in 
hard X-rays along the last five years. 

Our analysis also confirms that in the case of 
LS I +61°303 there is no high-energy cut-off found at ener- 
gies below 100 keV, as it would be common in conventional 
accretion scenarios. Indeed, if the system is an accreting neu- 
tron star or black hole, one expects to find a cut-off power- 
law spectrum in the hard X-ray band with a cut-off energy 
normally at 10 - 60 keV for neutron stars (e.g., Filippova 
et al. 2005) and at ~100 keV for black holes (McClintock 
& Remillard 2003). Even when one can think of situations 
where the accretion disk is masked by a jet component, pro- 
ducing a featureless power law, if the accretion flow and the 
jet would give comparable contributions in the hard X-ray 
band, one would still expect at least a feature in the IN- 
TEGRAL energy range, due to the disappearance of of con- 
tribution of the disk at high energies. Such is not observed, 
and thus, in case an accretion disk is present, its contribution 
should be sub-dominant with respect to the jet one. The flux 
we derive in the 20-200 keV band is 2.8 xl0 _11 ergs cm -2 



s~ . If most of it comes from a jet, and since in the scenario 
being discussion the whole source activity would be powered 
by accretion, the luminosity (Lj et ~ 10 35 erg s _1 , of the or- 
der of the gamma-ray luminosity) could only be explained 
by a large beaming factor. From ~50 mas resolution radio 
images of LS I +61°303 obtained with Multi-Element Ra- 
dio Linked Interferometer (MERLIN), extended, apparently 
precessing, radio emitting structures at angular extensions 
of 0.01 — 0.05 arcsec have been reported by Massi et al. (2001, 
2004). However, recent Very Long Baseline Array (VLBA) 
imaging obtained by Dhawan et al. (2006) over a full or- 
bit of LS I +61°303 has shown the radio emission to come 
from angular scales smaller than about 7 mas (which is a 
projected size of 14 AU at an assumed distance of 2 kpc). 
No large features or high- velocity flows were noted in any of 
the observing days, which implies at least a non-permanent 
jet nature. Changes within 3 hours were found to be in- 
significant, so the velocity of the outflow can not be much 
over 0.05c. The MAGIC collaboration also conducted a ra- 
dio campaign (in concurrence with TeV observations) to test 
these results using again MERLIN in the UK, the Euro- 
pean VLBI Network (EVN), and the VLBA in the USA 
(Albert et al. 2008). The results obtained by radio imaging 
at different angular scales show that the size of the radio 
emitting region of LS I +61°303 is constrained to be below 
~6 mas (~12 projected AU), and the presence of persistent 
jets above this scale is therefore excluded. As in the case of 
Dhawan et al. (2006), these observations have shown a radio- 
emitting region extending east-southeast from the brighter, 
unresolved emitting core. The outflow velocity implied by 
these observations is ~ 0.1c. These values of /3 would result 
too low to significantly increase the luminosity (Ljet oc S 4 ) 
even if strong jet precession is assumed. Finally, search for 
lines with the Chandra observatory yielded to no evidence 
for absorption or emission features in the source averaged 
and time-resolved spectra (Rea et al. 2010). Putting the IN- 
TEGRAL observations in this context, we find no obvious 
support for the accretion model in our data. 

As can be seen in Figure 4, the hard X-ray emission 
peaks in phase around 0.5, and in order to compare with the 
orbital phase derived at 7-rays, we have also generated the 



Long-term monitoring of LS I +61° 303 with INTEGRAL 5 



orbital light curves of 72XTE/PCA and .RXTE/ASM using 
the currently available observations. The orbital peaks of 
the soft and hard X-rays are close in phase (0.4-0.6) but not 
co-aligned. Thus, with low statistics, the hard X-rays may 
appear following the overall orbital modulation trend of the 
soft X-rays (e.g. see Harrison et al. 2000). We also show 
in Figure 4 the lightcurves obtained by Fermi (Abdo et al. 
2009, 8 months of observations since launch in June 2008), 
and MAGIC (several orbits put together, along few years 
of observations, as described in Albert et al. 2009). Direct 
inspection of these lightcurves seem to indicate that the hard 
X-ray component peaks in between the average soft X-ray 
and TeV lightcurves, but significant variability (e.g., orbit 
to orbit variability is already seen in MAGIC and VERITAS 
data, Albert et al. 2009, Acciari et al. 2009) may preclude 
to obtain conclusions in non-simultaneous observations. 
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